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Abstract 

We investigate in detail a recent model of colliding mobile agents [Phys. Rev. Lett. 96, 088702], used as an alternative approach 
to construct evolving networks of interactions formed by the collisions governed by suitable dynamical rules. The system of 
mobile agents evolves towards a quasi-stationary state which is, apart small fluctuations, well characterized by the density of 
the system and the residence time of the agents. The residence time defines a collision rate and by varying the collision rate, 
the system percolates at a critical value, with the emergence of a giant cluster whose critical exponents are the ones of two- 
dimensional percolation. Further, the degree and clustering coefficient distributions and the average path length show that the 
network associated with such a system presents non-trivial features which, depending on the collision rule, enables one not 
only to recover the main properties of standard networks, such as exponential, random and scale-free networks, but also to 
obtain other topological structures. Namely, we show a specific example where the obtained structure has topological features 
which characterize accurately the structure and evolution of social networks in different contexts, ranging from networks of 
acquaintances to networks of sexual contacts. 
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1. Introduction 

Till recently, graph theory has been used as the frame- 
work to study the dynamics and the topology of complex 
systems, where elementary constituents are represented 
by nodes and its interactions mapped to edges[l,2,3,4]. 
This approach enables the definition of statistical quan- 
tities, such as the clustering coefficient[5,6], the average 
path length and the degree distributions [3], which char- 
acterize the network structure and dynamics [2,4], be- 
ing applied to study physical [7], biological[8,9,10] and 
social[ll,12,13] networks. 

There are already many sets of real system data 
where these quantities were measured, e.g. data of 
virtual[14], social[15] and biological networks [16, 17], 
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and also data of large infrastructure networks[18,19]. 
However, while for some of these real situations graph 
theory describes the structure and the dynamical evo- 
lution of the underlying system[l,2,3], there are other 
situations where the standard approaches to construct 
networks do not agree with the underlying structure 
and dynamics [20,2 1,22]. A deeper understanding of 
these latter situations could be improved by consider- 
ing dynamical processes based on local information 
which produce such networks [6,11,23,24,25]. 

To this end, we proposed recently [24,25] a model 
where the nodes (agents) of the networks are mobile, 
with the connections being defined through their col- 
lisions, i.e. being a consequence of their dynamics. 
In other words, the network is build by keeping track 
of the collisions between agents, representing the in- 
teractions among them. Consequently, the network re- 
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suits directly from the time evolution of the system. 
With this approach we are able [24] to reproduce the es- 
sential features, namely the dynamical evolution, clus- 
tering and community structure, observed in empiri- 
cal networks of social contacts, and to explain associ- 
ated deviations of networks where the contacts are of 
a specific nature, e.g. sexual[24,25]. Social networks 
have attracted significant attention in the physics com- 
munity [11,12,13,26,27,28], where, using approaches 
and techniques from physical systems, gave further in- 
sight to understand particular social behavior and related 
phenomena, e.g. spread of epidemics [29,30], small- 
world effects [5,11], the mechanisms of social network 
growth [31] and cycles of acquaintances [6,32]. 

The main advantages of a model of mobile agents to 
model structure and evolution of social acquaintances 
comes from the fact that social interactions should be a 
result not of an a priori knowledge about the network 
structure but of local dynamics of the agents from which 
the complex networks emerge. Since those social inter- 
actions occur among individuals when they have close 
affinity or similar social features, it is understandable to 
consider networks of social contacts as the result of col- 
lisions among the agents moving in a continuous space 
or a projection of it whose associated metric is related 
to a social distance, defined from the sociologic factors 
promoting social contact. 

The aim of the present paper is twofold. On one hand 
it reviews the model in all its detail and discusses natu- 
ral ways to extend it for specific purposes. On the other 
hand it shows the features of the mobile agent model in 
what concerns quasi-stationary regimes, critical behav- 
ior and topological features. In Section 2 we describe 
the model of mobile agents. The model attains a quasi- 
stationary state which is fully characterized in Section 
3 and in Section 4 its topological and statistical features 
are presented. The system shows a percolation transi- 
tion which is described in Section 5. Application of the 
agent model to reproduce networks of social and sexual 
contacts is presented in Section 6 and finally Section 7 
discusses and concludes the main points throughout the 
text. 

2. Model of mobile agents with aging 

In this Section we will show that the model of mobile 
agents is parameterized by two single parameters, the 
density p characterizing its composition and the maxi- 
mal residence time controlling its evolution. As we 
will see, other quantities can be derived from these two 
parameters for particular purposes, for instance to study 




Fig. 1. (Color online) Snapshots of the system of mobile agents. 
Edges between two agents indicate that they already collided with 
each other: (a) Snapshot after the first collision and (b-c) two 
subsequent snapshots within one cluster of collisions. Filled nodes 
(blue) and unfilled nodes (red) represent two different types of nodes, 
e.g. males and females (see Sec. 6 for details). 

the critical behavior of the system. 

The system of mobile agents can be seen as a sort of 
granular gas[33], composed by N particles with small 
diameter d, representing agents randomly distributed in 
a two-dimensional system of linear size L ^ ^/N d and 
with low density p = N/ <C 1/d^. The mobile agents 
represent individuals moving through the system and 
each pair of agents may collide with each other, given 
rise to a social contact between them. This social contact 
is introduced by setting a link joining the two agents 
after the collision, till its removal when one of the agents 
leaves the system. Therefore, during the evolution of 
the system, each agent i is characterized by its number 
ki of Unks and by its age Ai . 

We consider periodic boundary conditions and when 
initialized each agent has a randomly chosen position 
and moving direction with velocity vq- When for the 
first time two agents collide, the corresponding collision 
is taken as the first contact in the network, as illustrated 
in Fig. la (see arrow). Letting the agents move and 
tracking their collisions through time, more and more 
connections appear (Figs, lb and Ic). Figure 1 shows 
the evolution of only one part of the system. Typically, 
several networks similar to the ones in Figs, lb and 
Ic appear during the evolution, each one formed from 
different single isolated collisions, as the one illustrated 
in Fig. la. 

The collision process is based on an event-driven al- 
gorithm, i.e. the simulation progresses by means of a 
time ordered sequence of collision events (diffusion) 
and between collisions each agent follows a ballistic 
trajectory [34] (drift). According to the algorithm, col- 
lisions take place whenever the distance between the 
centers of mass of two agents is equal to their diameter. 

Since collisions represent social contacts their dy- 
namical rules should fulfill some sociological require- 
ments. Namely, it is known[35] that many social inter- 
actions occur more commonly between individuals hav- 
ing akeady a large number of previous contacts. For in- 
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stance, in sexual contacts[25] individuals with a larger 
number of partners are more likely to get new partners. 
Therefore, we choose a collision rule where the veloc- 
ity of each agent may increase with the number k of 
previous contacts, namely 

where ki is the total number of social contacts of agent i, 
V = 1 m/s is a constant to assure dimensions of velocity, 
u) = cos 6 + By sin 6) with 9 a random angle and 
Ex and By are unit vectors. 

The exponent a in Eq. (1) controls the velocity up- 
date after each collision. For a = the velocity of each 
particle is constant in time, and consequently the kinetic 
energy density E = ^pv'^ of the system is constant. 
For a > the velocity increases with degree k. In this 
range, the value a = 1 {\v\ (x k) marks a transition be- 
tween a sub-linear regime (a < 1) and a supra-linear 
regime (a > 1) with different degree distributions[25]. 
Throughout the paper we will consider a = 1 in most 
of the situations, showing that it produces the suitable 
dynamics to reproduce real networks of social contacts. 
It should be noticed that while positive values of a yield 
dynamical laws which fulfill sociological requirements, 
the equation of motion (1) is also able to consider com- 
pletely different situations, where a < 0, i.e. where the 
ability to acquire new contacts decreases with the num- 
ber of previous contacts. In this manuscript, we will 
focus in the regimes for a > 0. 

As one may notice, contrary to collision inter- 
actions where the velocity vector is completely 
deterministic [36], here momentum is not conserved. 
This is a consequence not only of the increase of v 
but also of the fact that after one collision the moving 
direction is randomly selected. The main reason for 
this random choice is that, as a first approximation 
it is plausible to assume that social contacts do not 
determine which social contact will occur next. 

Concerning the residence time or 'age' parameter A, 
during which the agents remain in the system, if A ^ 
oo, each agent will eventually collide with all the other 
agents forming a fully connected network. Whereas, 
when the average residence time of the agents is fi- 
nite the system will reach a non-trivial quasi-stationary 
state[37,38], as described in the next Sec. 3. 

The aging scheme considered here is simply param- 
eterized by some threshold in the age of the agents. 
More precisely, each agent i is initialized with a cer- 
tain age Ai{0) which is a random number uniformly 
distributed in the interval [0, Tg] with Tg being the max- 
imal age an agent may have. Being updated according 
to Ai {t + At) = Ai (<) + At, the age eventually reaches 



Ai{t) = Tf, when the agent i leaves the system, yield- 
ing a total residence time Tf — ^^(0). Computationally 
the replacement of an old agent by a new one is carried 
out simply by removing all the connections of the old 
agent and updating its velocity to the initial value vq 
with a new random direction. 

When the average residence time is too small, two 
agents will have no time to collide at least once, and con- 
sequently no network is formed. On the contrary, when 
Te is too large, each agent will cross the entire system 
and a fully connected network appears. To avoid these 
two extreme regimes we consider an average residence 
time which is neither very small nor large when com- 
pared with the characteristic time r between collisions. 

For that, we define a collision rate, as the fraction 
between the average residence time Tf—(A(0)) = Tg/2, 
where {x) is the average of x over different snapshots 
in the quasi-stationary state, and the characteristic time 
T of the mean free path defined as 

\/2npd{v) ' 

where {v) is the average velocity of the agents. With 
this assumptions our collision rate reads 

A=^ = MZi. (3) 

where to is the characteristic time of the system at the 
beginning when all agents have velocity vq. 

Before proceeding in characterizing the behavior of 
such a system it is important to address three last points 
to understand the parallel between the model and real 
systems. First, the two-dimensional continuous space 
where nodes move is not the physical space where in- 
dividuals travel, meet or establish social acquaintances. 
Instead, it represents a projection of a highly dimen- 
sional Euclidean space whose metric is related to what 
is called the social distance[12,39]: the closer two nodes 
are, the similar their affinities are (same tastes, same be- 
havior, etc.) and therefore more probable to establish an 
acquaintance, i.e. to collide. It should be stressed that 
the metric is related to social distance, but also incorpo- 
rates effects of random factors promoting two persons 
to meet or establish friendship connections. We have no 
rigorous explanation for the fact that a two-dimensional 
projection of such a 'large' -dimensional space suffices 
to reproduce empirical data of acquaintances. But the 
fact is that it does, and therefore for simplicity we will 
consider a two-dimensional system. One could also con- 
sider higher dimensional systems of mobile agents but 
then the 'projected' velocity in Eq.(l) would change. 

Second, since the system of mobile agents is used to 
extract a complex network, the two parameters of the 
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model influence the statistical features of the network 
structure. Namely, increasing the density p confines the 
accessible region of agents thus promoting the occur- 
rence of collisions among them which are more con- 
fined in space. In other words, increasing the density 
one increases the clustering coefficient [5]. As for Ti, 
the larger the residence time, the larger the number of 
collisions an agent may have. Thus, increasing Tg in- 
creases the average number (fc) of connections. 

Finally, in a space of affinities what is the meaning 
of a velocity? The velocity is in fact a measure of the 
accessible region of a given agent within this social 
space which increases, decreases or remains constant 
after each collision, depending on the value of a. 



3. The quasi-stationary regime 

In this Section we describe in detail the non-trivial 
topological and dynamical properties of the quasi- 
stationary (QS) state of the model described in the 
previous section. As we will see, in the QS state the dy- 
namical and topological quantities fluctuate around an 
average value after a transient time of the order of 2Tt. 

Figure 2 illustrates the convergence toward the QS 
state for two different values of T^, namely for T( = 
30.75 and = 75.34. Here, the convergence is char- 
acterized by plotting the effective coordination given 
by the fraction M/N between the total number M of 
connections and the total number N of nodes (Fig. 2a), 
the mean energy {E) (Fig. 2b) and the mean age (A) 
(Fig. 2c), as function of time t. In all three cases, the 
above quantities increase in the earlier stages of the net- 
work growth attaining a maximum value around t ^ 
T(/2 where the agents start dying, resulting in the de- 
crease of their values till a minimum at i = Te. The 
maximum at t ^ Ti/ 2 is due to the fact that the agents 
have on average a life-time equal to /2, while the min- 
imum is due to the extinction of the first 'generation' of 
agents. Notice that, the maximum occurs at due to 
the choice of the initial uniform distribution of the node 
age which yields {A{0)); in general, for other choice of 
age distribution one finds a maximum at t ^ (yl(0)). 

The average values and corresponding fluctuations of 
these properties in the QS regime {t > 2Ti), depend on 
Ti, i.e. on A (see Eq. (3)). For large values of A (large 
Ti) the fluctuation is larger than for small A, as can 
be seen from the probability density functions (PDF) 
of each quantity shown on the right of Figs. 2a-c. The 
PDF of the data were computed from the first 2 x 10^ 
data points beyond t ^ 2Ti and are plotted with cir- 
cles, while a Gaussian with the same average and stan- 
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Fig. 2. Dynamical and topological quantities as function of time t: 
(a) effective coordination M/N, (b) mean energy (E), (c) average 
age (A). Two different values of the collision rate are plotted: 
Ti = 30.75 and Tf = 73.34. In (a) aiTows indicate t = Ti, the 
maximal age an agent may have. In all cases a quasi-stationary state 
is attained beyond t >2Ti. On the right the PDFs of these quantities 
are shown (circles), comparing them with Gaussians having the same 
central moments (solid lines). Here A'^ = 64 X 64 and a = 1 (see 
Eq. (D). 
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Fig. 3. Correlation fi between time step 2Ti, beyond which the 
QS is attained, and a subsequent time-step t > 2Ti, of the three 
quantities in Fig. 2: (a) effective coordination M/N, (b) mean 
energy {E), (c) average age (T). Solid line indicate the result for 
Tf = 30.75, while dashed lines indicate the result for = 73.34. 
Here iV = 64 X 64 and a = 1. 

dard deviation is plotted as a solid line for comparison. 
While for Ti = 30.75 the distributions are approximate 
Gaussian, for = 73.34 there is significant deviation 
from a Gaussian distribution. 

As explained below, between these two values there 
is a threshold corresponding to a critical collision rate 
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Fig. 4. Number G of nodes in the largest cluster of tlie system of 
mobile agents as a function of the collision rate A. At Ac ~ 2.04 
one sees the emergence of a giant cluster with G/N ^ 1 (see 
Sec. 5). Here Af = 128 X 128 and o = 1. 

Ac separating two different behaviors, one where the 
system is composed by several small clusters of inter- 
connected nodes and one where all the nodes belong 
to a single giant cluster. Further, from Fig. 3 one con- 
cludes that, beyond t > 2Ti, the correlation length is of 
the order of Tg. 

So, although there is no conservation of momentum, 
the average velocity increasing with time at the begin- 
ning, after attaining the QS state {v) is almost constant. 
In the following we assume only values of A in the QS 
state, taking the average velocity in Eq. (2) as a constant 
given by the mean of the corresponding PDF distribu- 
tion. Further, although all quantities, M/N, (E) and 
(T), characterizing the system in the QS state are func- 
tions of Te, we will choose the adimensional parameter 
T(/to instead, which incorporates also the density of 
the system. 

As said above, by varying the collision rate A, one 
finds a critical value Ac marking a transition from a state 
composed by several small clusters to a state where a 
giant cluster emerges after attaining the QS state. In 
Fig. 4 we plot the fraction G/N between the number 
G of nodes in the largest cluster and the total number 
N of nodes. Clearly, beyond Ac — 2.04 one sees the 
emergence of a giant cluster. In the next Section 4 we 
describe and discuss the topological properties for this 
largest cluster and in Section 5 we characterize the per- 
colation transition occurring at Ac. 

4. Properties of the network 

In this section we will study the degree distribution 
P{k), the clustering distributions C(fc) characterizing 
the cliquishness of the agents neighborhoods, and the 
average path length i in the QS state. 




Fig. 5. Average degree (fc) as a function of (a) the collision rate A 
and of (b) the maximal age T^/tq for two values of the density, 
p = 0.02 (white circles) and p = 0.2 (black circles), both with 
a = 1. The solid line in (a) indicates A = 2{fc). (c) Degree 
distribution of the giant cluster in the QS state for two different 
velocity updates (see Eq. (1)): a = (triangles) with Ti/tq = 9.5, 
where the velocity modulus is always constant and o = 1 (diamonds) 
with Tf/ro = 3, where velocity modulus increases hnearly with 
the number of previous connections. Measuring in each case (fc) in 
the system and introducing its value in Eqs. (4) and (5), yields the 
suitable fits for each a value, namely a Poisson distribution (solid 
line) with (fc) = 2.52 and an exponential distribution (dashed line) 
with (fc) = 2.62. In all cases, averages over 100 iterations were 
taken after attaining the QS state for a system of N = 10** agents. 

The average degree (k) per agent is a function of both 
the collision rate A and rescaled maximal age T^/tq as 
shown for two different values of p in Fig. 5a and 5b 
respectively. From Fig. 5a one sees that (fc) = A/2 in- 
dependently on the density p, while as a function of 
Ti/tq the average degree increases exponentially with 
an exponent depending on the density. These two plots 
are important not only to study the topological features 
of the agent model, but also to find the appropriate val- 
ues when modeUng real networks as will be explained 
in Section 6. 

In Fig. 5c we compute the degree distribution P{k), 
by counting the fraction of nodes having k neighbors. 
Two different velocity updates are illustrated here, 
namely a = (triangles) and a = 1 (diamonds). 
Clearly, the degree distribution depends strongly on the 
collision update rule, i.e. on the value of a in Eq. (1). 

More precisely, for a = 0, the velocity of each agent 
is always constant, the resulting degree distribution be- 
ing a consequence of the effective time [24] t^ff = 
3T( /8 to create links and of the collision rate, yielding 
a Poisson distribution 

PW = ^e-<^>. (4) 

For this value a = 0, the degree distribution was cal- 
culated for a fixed value of T;/to = 9.5 the resulting 



5 




1000 



Fig. 6. (a) Degree distribution of tlie giant cluster in the QS state, for 
several values of A = 5.24, 8.04, 13.04, 23.76 and 57.36 (symbols). 
Lines indicate the con'esponding exponential tit with Eq. (5). Here 
A'^ = 64 X 64 and a = 1. For other system sizes the results are 
similar. 
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Fig. 7. Clustering coefficient distribution for the same simulations 
of Fig. 6, namely for (left to right) A = 5.24,8.04, 13.04,23.76 
and 57.36. The solid line is a guide to the eye with slope 1/2 (see 
details in text). Here a = I, N = 128 X 128 and averages over 
100 realizations were taken. 

network has (fc) = 2.52, introducing this value in eq.4, 
we obtain the solid line in Fig. 5c, which shows that 
the degree distribution of the network is well aproxi- 
mated by a poissonian. In this way, one can argue that 
for a = the system produces a two-dimensional geo- 
metric random graph [40] in the QS state. 

For a ~ 1 the velocity in Eq. (1) increases linearly 
with the number of previous collisions. As we will see, 
this kind of dynamics reveals to be most suited to repro- 
duce the statistical features of real social networks (see 
Sec. 6). In this case, the effective residence time is uni- 
formly distributed (see Sec. 2), yielding an exponential 



velocity distribution and consequently an exponential 
degree distribution 

1 



P(fc) = 



(5) 



(fc)-l 

In Fig. 5 the dashed line indicates the distribution in 
Eq. (5) with (k) = 2.62, which results from a value of 
Ti/tq = 3.0 used in the numerical simulation (trian- 
gles). So, for this value a = 1 the agent model is capa- 
ble of producing a two-dimensional exponential graph 
in the QS state. 

In the following we will consider the value a = 1 in 
Eq. (1) and study how the degree distribution depends 
on the collision rate A. Figure 6 shows the degree dis- 
tribution in the QS state for A = 5.24 (circles), 8.04 
(squares), 13.04 (diamonds), 23.76 (up triangles) and 
57.36 (left triangles), plotting with lines the exponential 
in Eq. (5) evaluated with the corresponding (fc). As one 
sees, while for small A the exponential expression fits 
well the observed degree distributions, for larger (fc) the 
numerical results have a lower cutoff than the analyti- 
cal expression. Therefore, one concludes that the agent 
model is able to reproduce exponential distributions for 
low values of the collision rate (A < 10), and that other 
non-trivial distributions appear increasing the collision 
rate. The latest ones are the ones observed in empirical 
social networks (see Sec. 6). 

Another property of interest to characterize the net- 
work is its clustering coefficient. The local clustering 
coefficient C; of a vertex i with degree fc; measures 
its cliquishness and is defined as the quotient between 
the number Wi of triangles (cycles composed by three 
edges) and the total number of triangles. 

To uncover hierarchical properties of the network, one 
usually [41] computes the clustering coefficient as a 
function of the degree fc, namely 

2{wik)) 



C{k) 



(7) 



fc(fc- 1)' 

where {w{k)) is the average number of triangles of a 
node with degree fc. 

Figure 7 shows the mean degree-dependent clustering 
coefficient C(fc) for the same mobile agents systems of 
Fig. 6. It is interesting to note that, in contrast to random 
graphs, which have a clustering coefficient independent 
of fc, here we observe a dependence of the form ^ fc~" 
with a G [0.4,0.6]. In other words. Fig. (8) shows that 
the clustering coefficient decreases with the degree, a 
feature which indicates the existence of an underlying 
hierarchical structure [42]. 
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Fig. 8. Average path length £ (•) and clustering coefficients C (o) 
as functions of the average degree. Here p = 0.02. 

To end this Section, we study two other topological 
quantities, namely the average shortest path length £ and 
the average clustering coefficient C. As one knows [5], 
knowing the values of these topological quantities one 
is able to ascertain for small-world effects. 

The average clustering coefficient is defined as the 
average of the local clustering coefficients Ci in Eq. (6) 
for alH = 1, . . . , TV, and the average shortest distance £ 
is given by the average minimal number of edges joining 
two randomly chosen nodes. 

To investigate small-world effects in our system, we 
compare the networks of the model of mobile agents 
with random networks having the same number N of 
nodes and where each pair of nodes is connected with a 
probability p = 2M/{N{N~1)). With this probabiHty 
one obtains a random network with approximately the 
same effective coordination M/N (same (fc)) as the one 
observed in the model of mobile agents. 

Figure 8 shows these two quantities as a function of 
the average degree (fc). The clustering coefficient in- 
creases with (k), and beyond the emergence of the giant 
component ((fc) > Ac/2 = kc) it slowly converges to 
a value Coo ^ 0.1. The vertical line shows the critical 
point kc = 1.02. As for the average path length £ one 
observes a slow decrease with the average degree. No- 
tice that for very low values of (fc) (not shown) there 
are no sufficient edges to guarantee a meaningful aver- 
age path length, in particular, below the transition at kc 
the system is divided in several different clusters, which 
leads to an undefined £. 

In Fig. 9a, the average path length is small compared 
to the system size. Moreover, as seen in Fig. 9b the 
cluster coefficient in the agent model (circles) is much 
larger than in the random counterparts (triangles). How- 
ever, the networks generated by the model of mobile 
agents are not small-world, since to have the small- 
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Fig. 9. Size dependence of the networks of mobile agents, (a) 
Shortest path length ^ as a function of the system size N, for 
Ti/to = 5.02, compared to In A'^ (dashed lines). For large network 
sizes the sohd line indicates a fit ^ = e~^V~N (see text), (b) 
Average clustering coefficient (C) as a function of A'^ for the agent 
model (bullets) compared to the corresponding random graph with 
the same (fc) (triangles), having a fit (solid line) of C ~ 1/N. Here 
p = 0.02. 

world property[3] it is also required that the increase of 
the shortest path length with the system size is not faster 
than IniV. From Fig. 9a one sees that this is only true 
for small system sizes {N < 10"*). For larger systems, 
the fitted numerical results yield £ ~ e^'^^fN. 

In Fig. 9b one sees the behavior of the average cluster 
coefficient when the system size is increased. Interest- 
ingly, one clearly sees an independence on N, beyond 
N > 10'^, where C ~ 0.08. For higher density values 
p this A^-independent value of C increases. This result 
is quite in contrast to random graphs which scales with 
1/N, as illustrated in Fig. 9b with triangles. 

Having characterized the features of the system be- 
yond the transition at Ac, where a giant cluster emerges, 
in the next Section we focus on the properties of this 
transition. 

5. Critical behavior 

The emergence of the giant cluster for an exponential 
network has been studied by Molloy and Reed [43] or 
equivalently using the generating function formalism 
by Newman [44] . In both cases the transition belongs 
to the universality class of mean field percolation, in 
contrast to the agents model which is in the universality 
class of two-dimensional percolation, as we will show 
in this Section. 

The percolation transition is the abrupt change, at a 
particular critical value Ac, from a state composed by 
several small clusters to a state where a giant cluster 
dominates. A cluster is defined as a group of connected 
agents. Isolated agents are regarded as clusters of size 
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unity and any cluster consisting of s connected agents is 
called an s-cluster. We borrow the notation from Stauf- 
fer's book on percolation theory [45] and define here 
71 s = Ns/N as the number of s-clusters per agent, 
where Ns is the number of clusters of size s and N the 
total number of agents in the system. 

In Fig. 10 we show the time evolution of the size G 
of the largest cluster (see also Fig. 4), together with the 
total number ^?is of clusters and the corresponding 
size-averaged mean cluster size ^ s^Tig. At the begin- 
ning the size of the giant cluster and the size of the 
mean cluster increase till t = T(/2, corresponding to 
a complementary decrease of the number of clusters. 
Beyond this instant, agents which have an average res- 
idence time of start to leave the system losing all 
their connections, and therefore both the size of the gi- 
ant cluster and the mean cluster size start to decrease 
till t ^ Tg. Thereafter, all the first generation of nodes 
has left the system, and the evolution converges rapidly 
to the QS state where all quantities present constant 
values, apart small fluctuations. For other values of T^, 
i.e. of the collision rate A, one observes the same be- 
havior of the above quantities, reaching the QS state at 
different values. 

The next step is to calculate for a given value of T;, 
the value of A in the QS state and average at different 
times the value of G. The result is plotted in Fig. 11a 
for different system sizes. Simultaneously, the value of 
X, given by the size-averaged mean cluster size without 
the largest cluster, is also calculated. The result is shown 
in Fig. 12a. 

In general, according to the standard scaling the- 
ory [46], we expect G (Fig. 11a) and x (Fig- 12a) to 
follow 
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Fig. 10. Cluster quantities as function of time t: Number G of 
agents belonging to the largest cluster (solid line), total number 
of clusters SsTIs (dotted line) and size-averaged mean cluster size 
TisS^Us (dashed line). All quantities are normalized by A'^. Here 
AT = 128 X 128 and A = 13.04. 




A, (X-XJL'" 

Fig. 11. (a) Fractions of sites in the largest cluster G/N as function 
of A. The numerical results where extracted in the QS state for 
a fixed value of Tj/tq, with N = 322,46^,922, 128^ and 256^. 
(b) Confirmation of the scaling relation in Eq. (8), for the system 
sizes A'^ = 2^", 2^^, . . . , 2^'' and the values of /3 and v reported 
in Tab. 1 at the critical value Ac = 2.04. 




Fig. 12. (a) Size-averaged mean cluster size ^ s^ns as a function 
of A for the same system sizes as in Fig. 11a. (b) Confirmation of 
the scaling relation of Eq. (9), for the same system sizes the values 
of 7 and u reported in Tab. 1 and Ac = 2.04. 

G' = L-f/''i^[(A- AjL-i/"^] (8) 

X = L''/''G[{\-X,)L-^/n- (9) 
where L is the linear size of the system {N = L^). 
This is confirmed by the collapses of the curves near 
the the critical point, i.e. A — Af, ~ 0, as illustrated in 
Figs, lib and 12b respectively, for the values of i^, j3, 
and 7 reported in the central column of Tab. 1, with 
Ac = 2.04. 

In order to test the scaling relation obtained above, 
we compute the number ns(A) of clusters of size s as 
a function of A, as shown in the inset of Fig. 13, for 

= 2^^ agents. For the computation, we calculate bins 
(power-two) of the cluster size distribution at different 
values of A, ignoring the first four bins, that is the sizes. 
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MF mobile agents Percolation (2D)[45] 
u 0.5 1.3 ±0.1 4/3 ~ 1.33 
7 1 2.4 ±0.1 43/18 ~ 2.39 
(3 1 0.13 ±0.01 5/36 ~ 0.139 
0-0.5 0.40 ±0.01 36/91 ~ 0.397 

Table 1 

Critical exponents related to the emergence of the giant cluster in a 
random graph model (percolation MF), for the network of mobile 
agents presented here, compared to the exact results of percolation 
2D. 

1, 2 — 3, 4 — 7, and 8 — 15, since for such small clusters 
scaling is not good[45]. We take the next 4 bins, i.e. s- 
values in the ranges [16 - 31], [32 - 63], [64 - 127] and 
[128 - 255] and plot them for different A. 

The main plot in Fig. 13 shows the result after scaling. 
Here, for each bin we take the ratio ns{\)/ns{\c) as a 
function of (A — Ac).s'^ using s as the geometric average 
over the two extremes of the bin ranges enumerated 
above. Performing this scaling we obtain 



1 



/3 + 7' 



(10) 



using the values of jS, 7 and Ac obtained previously (see 
Tab. 1). 

Equation (10), confirms that the scaling relation of 
the phase transition belongs to the universality class of 
percolation. It is known that the emergence of a giant 
cluster for a random graph depends on (fc) with its crit- 
ical value at (fc)^ = 1, and the phase transition belongs 
in the same universality class as mean field percola- 
tion [47], whose exponents are given in the first col- 
umn of Tab. 1 . However, for the agent model one ob- 
serves the same emergence of the giant cluster at (fc)c = 
0.5 Ac = 1.02 and the universality class corresponds to 
2(i-percolation, as shown in Tab. 1 . 

Before ending this Section, we stress that the corre- 
lation exponent v presented in Tab. 1 is obtained from 
finite size scaling. This exponent can be also explicitly 
calculated from the linear size of clusters (see [45]), 
namely by computing the correlation length ^(A) ^ 
|A — Acl""^. The result is shown in Fig. 14, where the 
solid line has a slope of —1.3 yielding a correlation 
exponent in agreement with the previous results (see 
Tab. 1). Since the agents move on a two-dimensional 
plane and have only a finite life time, they can only 
establish connections within a restricted vicinity, and 
this effect corresponds to a connectivity which is short 
range. Thus, although the clusters in the agent model are 
not quenched in time, the underlying dynamics yields a 
short range 2(i-percolation. 




Fig. 13. Collapse of the curves from the inset, plotting ns/ns{\a) 
as a function of (A — \c)s" (see text). In both cases curves were 
plotted for cluster sizes in the ranges [16,31], [32-63], [64-127], 
and [128 - 255]. 
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Fig. 14. CoiTelation length § as function of A — Ac. Symbols indicate 
the result of simulations performed for different values of A and the 
solid line has a slope of v = —1.3. 

6. Real network of social interactions 



In the previous Section we described the main char- 
acteristics of the model of mobile agents in the QS 
regime, namely an exponential degree distribution and 
a transition to percolation for a certain critical collision 
rate, beyond which no small-world effects are observed. 
In this Section we analyze a real social network and 
show that its topological features are well reproduced 
by our model. In particular, we show that the statistical 
features observed for the empirical networks are well 
reproduced with parameter values beyond the transition 
to percolation. 

The empirical data comprehend an extensive study 
done within the National Longitudinal Study of Ado- 
lescent Health (AddHealth) [48] at the Carolina Popula- 
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Fig. 15. Reproducing tlie distribution of friendsliip acquaintances in 
empirical social networks [48] with the agent model. In all cases 
circles represent the empirical data, while solid hues indicate the 
result of the model of mobile agents. 

tion Center, concerning 84 American schools. The aim 
is to allow social network researchers interested in gen- 
eral structural properties of friendship networks to study 
the structural and topological properties of social net- 
works [32]. All data are constructed from the in-school 
questionnaire from a total of 90118 students which re- 
sponded to it in a survey between 1994 and 1995. The 
students are separated by the school they belong to in a 
total of 84 schools whose corresponding network sizes 
range from ~ 100 to ~ 2000 students. 

Figure 15 shows several plots of the degree distribu- 
tions in nine of the schools (circles) compared to the 
distributions obtained with the agent model (solid lines) 
for the same number of agents. The values for the col- 
lision rate in the simulations were taken from A = 2(fc) 
where the average degree is the one found in the cor- 
responding school. As one sees, in all cases the tails of 
the distributions are well fitted by the model. 

Clearly, the agent model reproduces the degree dis- 
tribution observed in schools, by introducing the same 
average degree found in real data. The agent model is 
not only able to reproduce these degree distributions, 
shown in Fig. 15, but also the corresponding clustering 
coefficient and the shortest path length. 

Figure 16a shows for each schools (symbols) the av- 
erage shortest path length £ (circles), and the clustering 
coefficient (triangles). Solid lines indicate the results 
obtained for the agent model using the same values of 
(fc), averaged over 100 realizations with N = 2209 
and p = 0.1. Since d depends on the network size, it 
is divided by the shortest path length Iq of a random 
graph with the same average degree and size. Clearly, 
the agent model predicts accurately both the clustering 
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Fig. 16. (Color online) (a) Average shortest path length I and clus- 
tering coefficient C as functions of the average degree (fe). Empiri- 
cal data (symbols) compared to simulations (solid hues), (b) Plot of 
Til To as a function of (fc) for the agents models (solid fine). Stars 
illustrate two particular schools for Figs. 15 having Til To = 4.75 
(school 1) and 6.0 (school 2) respectively, (c) Second moment (fc^) 
for each school vs. the second moment of the corresponding simu- 
lation with the agent model (solid fine has slope one). In all cases 
p = 0.1. 

coefficient and the shortest path length for the same av- 
erage degree. 

By computing the average degree {k) of each school 
one is able to obtain the value of T/ / tq for which the 
agent model reproduces properly the empirical data, as 
illustrated in Fig. 16b. Here solid lines indicate the cal- 
culated curve from the agent model, while triangles in- 
dicate the values of Ti /tq chosen to reproduce the social 
network of the schools with the resulting value of {k). 
Moreover, the second moment {k'^)ag obtained with the 
simulations of the agent model is rescaled by the same 
quantity {k'^)sch measured for the empirical school net- 
works, as shown in Fig. 16c. 

In other words, a strong point of our model is that 
since the average degree varies monotonically with A 
one can associate a particular value of the collision rate 
to each school. The collision rate and in particular this 
parameterization for the collision rate depends weakly 
on the network size. Here we neglect this dependence, as 
a first approximation. Using this one-to-one correspon- 
dence one is able to compare the topological quantities 
between the schools and the model of mobile agents. 

Notice that, the Add Health in-school friendship nom- 
ination data was constructed from a questionnaire where 
each student listed their in-school friends by order of 
importance (best friend first) and to each of the friends 
they indicate which kinds of contacts they had in the 
last seven days, out from five possible choices. Accord- 
ing to the number of different kinds of acquaintances a 
weight from 1 to 5 is attributed to each social contact. 
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Fig. 17. Cumulative degree distribution of the number k of sexual 
partners in a real empirical network of sexual contacts (triangles) 
with 250 individuals, compared to the simulation of the agent model 
(solid line), the dotted line is a guide to the eye with slope 2. Here 
N = 4096, Ti/to = 5.5 and (fc) = 7.32 and the average size of 
the resulting sexual network is 220. 

Moreover, since the questionnaires are based on the in- 
formation given by one agent, the networks are directed. 
Here, we each social contact. Moreover we assume as a 
first approximation that social contacts are unweighted. 
Using all the information concerning these empirical 
networks [32] a more detailed study of the structure of 
these networks will be presented elsewhere [49]. 

Finally, we stress that with a shght extension of the 
model it is possible to reproduce subsets of social con- 
nections, e.g. sexual contacts, having particular features 
different from the entire network of all the acquain- 
tances. For instance, in the case of sexual contacts the 
degree distribution is commonly a power-law [21]. Fig- 
ure 17 shows the degree distribution of a sexual contact 
network (triangles) extracted from a tracing study for 
HIV tests in Colorado Springs [50]. The dashed line in- 
dicates the result obtained for a network of all social 
contacts simulated with the agent model, while the solid 
line is the degree distribution of a subset of contacts 
from the social network. 

The contacts in the subset are chosen by assigning 
to each agent an intrinsic property, say a fitness, which 
enables one to select from all the social contacts the 
ones of particular interest, in this case sexual contacts. 
The fitness is defined by a given number and, when 
two agents form a link, this link is marked as a 'sexual 
contact' if the sum of agents fitnesses, is greater 

than a given threshold. The initial distribution of the 
fitness is assigned to the agents following an exponential 
distribution and the threshold is lnA^/2, following the 
scheme of intrinsic fitness proposed in another context 
by Caldarelli et. al. [51]. Interestingly, one is able to 
extract from the typical distributions of social contacts 
shown in this Section, which are not scale-free, power- 
law distributions in QS state which resemble much the 



ones observed in real networks of sexual contacts (see 
also Ref. [25]). 

7. Conclusions 

In this paper we introduce an approach to construct 
networks of complex interactions from a model of mo- 
bile agents. The network is constructed by keeping track 
of the collision between the mobile agents. Introducing 
an aging scheme the network attains a QS regime, char- 
acterized by almost constant topological and dynami- 
cal quantities which fluctuate slightly around an aver- 
age value. Typically, the networks constructed with the 
agent model, while having small path lengths and large 
clustering coefficients, do not show small-world effects, 
since the average path length does not scale logarithmi- 
cally with N. 

The control parameters of the model are the density 
of the agents and the maximal residence time the agents 
may have. With these two parameters it is possible to 
define a collision rate which, for a critical value, marks 
a percolation transition whose exponents are the ones 
of two-dimensional percolation. While for small colli- 
sion rates, the degree distributions are approximately 
exponential, for larger values of the collision rate they 
strongly deviate from the exponential approximation 
yielding a non-trivial degree distribution. Surprisingly, 
these degree distributions are precisely the same as the 
ones observed in real social networks of empirical data. 
This empirical data has many other informations which 
were not used in this study, namely the gender of the 
agent and the 'direction' and 'weight' of the acquain- 
tance. Therefore future investigations on the context of 
network analysis of this recent empirical data could be 
interesting for sociological purposes [49]. 

Moreover, we have shown that the agent model is also 
able to reproduce networks having only social contacts 
of particular nature, e.g. sexual contacts. In the partic- 
ular case of sexual contacts it was found [21] that the 
distribution of sexual partners has a power-law tail. The 
dynamical reasons underlying this discrepancy between 
the 'total' number of social contacts, with exponential 
distributions, and the particular number of sexual con- 
tacts, with power-law tails, is not yet explained. 
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